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The definition and measurement of magnetic reconnection in tliree-dimensional magnetic fields witli multiple 
reconnection sites is a challenging problem, particularly in fields lacking null points. We propose a general- 
ization of the familiar two-dimensional concept of a magnetic fiux function to the case of a three-dimensional 
field connecting two planar boundaries. Using hyperbolic fixed points of the field line mapping, and their 
global stable and unstable manifolds, we define a unique flux partition of the magnetic field. This partition 
is more complicated than the corresponding (well-known) construction in a two-dimensional field, owing to 
the possibility of heteroclinic points and chaotic magnetic regions. Nevertheless, we show how the partition 
reconnection rate is readily measured with the generalized flux function. We relate our partition reconnection 
rate to the common deflnition of three-dimensional reconnection in terms of integrated parallel electric fleld. 
An analytical example demonstrates the theory, and shows how the flux partition responds to an isolated 
reconnection event. 



I. INTRODUCTION 

This paper presents a new method for measuring mag- 
netic reconnection in a three-dimensional (3D) magnetic 
field. Reconnection is a fundamental physical process in 
any highly-conducting plasma, yet remains poorly un- 
derstood owing to the challenging range of lengthscales 
involvedii^. In 3D magnetic fields, progress is hampered 
by the difficulty in defining and measuring reconnected 
fiux, particularly if there are multiple interacting recon- 
nection sites. It is this problem that we seek to address. 

There are two contrasting ways to measure reconnec- 
tion rates in 3D. The first uses the parallel electric field 
integrated along magnetic field lines^, while the second 
counts the transfer of fiux between distinct flux domains. 
Here we pursue the second approach, where the task is 
twofold: to define a partition of the flux and to measuring 
the rate of transfer between fluxes in this partition. We 
call this a reconnection rate with respect to a partition, 
or a partition reconnection rate, to distinguish it from 
the first case. Such a partition reconnection rate can 
capture only reconnection processes that change fluxes 
between the partition domains, and not those within any 
individual flux domain. However, this rate is in many 
applications the most relevant information, detemining 
the stability and dynamics of the system. 

In a two-dimensional (2D) field B = Bx{x,y)e.j; + 
By{x,y)ey, there is a natural choice for such a partition 
and correspondingly for the reconnection rate: write B 
in terms of a flux function A{x, y) where B = V x Aez- 
The different fluxes in the partition correspond to the re- 
gions of the plane bounded by separatrices, which are the 
global stable and unstable manifolds of hyperbolic nulls 
(x-points). The magnetic flux within each such region 
(per unit height in the ignorable z direction) is measured 
by the difference in A between appropriately chosen nulls. 
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These fluxes are invariant under an ideal evolution, while 
in a non-ideal evolution the change in fluxes is measured 
exactly by the change in the values of A at the (discrete) 
set of null points. This deflnes an unambiguous global 
reconnection rate which is readily computed even in tur- 
bulent 2D fields with many nulls^. 

In 3D, the situation is more complicated. Here a natu- 
ral partition also arises from the existence of null points 
in the domain. The 2D invariant manifolds (fan sur- 
faces) associated with these null points form a coarse 
but natural partition of the fluj*^!^. This inherent topo- 
logical structure has been used successfully to quantify 
reconnection^. Nulls and separators in particular are 
favoured locations to detect global changes in connectiv- 
ity because they are locations where distinct flux domains 
come into close proximity. 

Another possible flux partition arises in a fleld con- 
nected to a physical boundary, such as the photosphere 
of the Sun, where the sign of the normal field component 
divides the boundary into regions of positive and nega- 
tive magnetic polarity. This boundary partition extends 
to a partition of magnetic flux connected to the bound- 
ary, by following the field lines from the different polarity 
regions into the volume^. 

There are, however, many examples of magnetic fields 
where either the above described partitions are too coarse 
or null points do not exist in the domain. Examples are a 
single coronal loop or the magnetic field in a tokamak. A 
generic 3D magnetic field in this situation does not pos- 
sess a foliation of flux surfaces. This "worst case" is the 
situation which we want to investigate here. More specif- 
ically, we assume a simply-connected domain in which 
all fleld lines stretch between two planar boundaries. In 
magnetospheric reconnection studies this is often referred 
to as the "guide field" case. We present a general way 
to define the flux partition in such a fleld, using distin- 
guished hyperbolic orbits, and measure reconnection by 
introducing a generalized version of the 2D flux function 
A. Not only is this a natural topological flux partition 
when there are no magnetic null points, but it retains the 
simplicity of the 2D method when it comes to measuring 



the partition reconnection rate. 

As may be expected, there are some complexities that 
do not arise in the 2D case. Chief among these is 
the possibihty of chaos in the field line mapping. This 
is well-known from the study of area-preserving map- 
pings as models for the magnetic field in toroidal fu- 
sion devices^. Indeed, recent results in tokamak exper- 
iments show heating on the vessel walls consistent with 
the breakdown of confinement and chaotic transport of 
magnetic flux through homoclinic tangles, as found in nu- 
merical siniulationsiSrJ^. Though they were recognised 
by Poincare in the 19th Century, it is only recently that 
detailed analysis of the structure of homoclinic tangles 
has been applied to measure and predict the transport 
of trajectories in these chaotic regions. A primary appli- 
cation has been 2D fluids with time-dependent velocity 
field^i^"— . Here, we show how these ideas can be applied 
to define and measure a natural reconnection rate in 3D 
magnetic fields. 



II. TWO-DIMENSIONAL MAGNETIC FIELDS 

We first review the basic properties of the flux function 
A{x,y) of a 2D magnetic field B = V x Ae^ (Figure [1]). 

1. A is constant along magnetic field lines (B • VA = 
0). 

2. Consider two vertical lines through the points 
{xi,yi) and {x2,y2)- The magnetic fiux through 
any surface bounded by the two lines and the planes 
z = 0, z = 1 is A(x2-, j/2) — A(xi, j/i) (w.r.t. the ori- 
entation of the surface). This is the reason for the 
name flux function. It works because 



A{,x,y)^ / 
Jo 



A{x, y)dz = / A- d\ 



(1) 



where A = A{x, y)ez is a vector potential for B. 
3. For an ideal evolution 
9B(x,y,i) 



dt 



Vx (v(a-,j/,t)xB(x,y,i)) == 0, 



A(x,y) can be chosen as an ideal invariant, 

^+v.V^^O. (2) 

A 2D magnetic field is naturally partitioned by the 
x-points (hyperbolic nulls) and their separatrices. The 
scparatrices — shown by thick lines in Figure [1] — are the 
topologically distinguished field lines given by the global 
stable and unstable manifolds of each x-point. These 
manifolds are tangent at the null to the unstable or stable 
eigenvectors of the local linearisation, and are uniquely 
defined by extending forwards or backwards along the 
fiow. They are invariant subspaccs of the field line flow 



(i.e., they are field lines) that delineate topologically dis- 
tinct regions. Each region has a well-defined flux mea- 
sured by the difference in A between two (not necessarily 
unique) null points (joined by dashed lines in Figure [!}. 
Since the separatrices are themselves field lines, two x- 
points joined by a separatrix must have the same value, 
of A. 




FIG. 1. A 2D magnetic field, showing nulls (hyperbolic 
as squares, elliptic as circles) and separatrices (thick lines). 
Dashed grey lines show the differences in A that measure the 
partition fluxes. 



In a 2D field, changes in topology — i.e., changes in 
the amount of fiux in each region of the partition — can 
take place only at null points^i. Traditionally, recon- 
nection counts only changes in A at x-points, indicat- 
ing the transfer of fiux between distinct regions, and not 
at o-points, where changes in A represent only the cre- 
ation/annihilation of flux within a single region. In this 
way one can define a global partition reconnection rate 



A$t: 



E 



d^(h,) 



dt 



(3) 



where the sum is over all x-points h^. Notice that we can 
measure the reconnection rate completely knowing only 
the values of A at the null points, with no need to know 
either the partition fiuxes or the structure of the field 
(except for the spatial derivatives at each null point, in 
order to determine the hyperbolicity) . 



III. GENERALIZED FLUX FUNCTION 

In a general 3D field, we can no longer write B in 
terms of a 2D function. But this does not prevent us 
from constructing a 2D function to measure magnetic 
flux. In this section, we will show how to construct a 
generalized flux function A{x^ y) that retains a number 
of the properties of the 2D flux function Aix^y). Our 
domain is a bounded region in 'W? between z = and 
z = 1, with all field lines connecting from the lower to 
the upper boundary. 



As a simple generalization of A{x,y), we might con- 
sider a function f{x,y) = J^ A ■ ezdz, i.e. the inte- 
gral along vertical lines of the vector potential A (where 
B = V X A) . The difference in f{x, y) between two points 
(xi,yi) and (a;2,J/2) would then give the flux through a 
vertical surface, analogous to the 2D case. However, this 
function f(x,y) does not retain the ideal invariant prop- 
erty of A{x,y), which is vital to define any meaningful 
reconnection rate. 

To construct a flux function that is an ideal invariant, 
we make a simple modification and integrate A along 
magnetic field lines rather than vertical lines. For a point 
(a;, y) on the lower boundary, denote the field line starting 
at {x,y) by Fz{x,y). In other words. 



dF,{x,y) B(F,(x,j/)) 



with Fo{x,y) ^{x,y). 



dz B4F4x,y))' 

(4) 
The subscript z indicates that we have chosen to 
parametrise the field line by the vertical coordinate z. 
With this notation, we may define the generalized flux 
function as a function on the lower boundary z = 0: 



A{x,y) 



A(F.(.,„))-|§Hr''- 

B,{F,{x,y)) 



(5) 



We consider in this paper only periodic fields where 
Bz{x,y,l) — Bz{x,y,0), and we impose the gauge con- 
dition that A X e^ is periodic. We are still free to impose 
a gauge transformation A — > A' -I- Vx providing that 
x{x, y, 1) = x(a;, y, 0) -|- xo with xo constant. We impose 
the further gauge condition xo = 0; leaving the function 
xi^i 2/i z) free for < z < 1. 

Under a gauge transformation, the function A{x, y) be- 
comes 

A'{x,y) ^ A{x,y) + x{Fiix,y)) - xix,y,0), (6) 

so it is not gauge invariant in general. But at fixed points, 
where Fi{x,y) = {x,y), the last two terms in ^ cancel 
and A{x, y) becomes gauge invariant. Thus diff'erences in 
A between fixed points are well-defined, and correspond 
to physical fluxes (Figure [2]), in analogy to the 2D case. 
We argue in this paper that these physical fluxes defined 
by values of A at fixed points form a natural partition of 
the 3D magnetic field. In fact, the values of A at non- 
fixed points may also be given physical meaning if one 
fixes the gauge in a particular way; this is beyond the 
scope of this paper and will be addressed in future. 




FIG. 2. The flux $ioop through the surface deflned by two 
fixed points {xi,yi), (2:2,2/2) is measured by A{x2,y2) — 
A{xi,yi), because the integrals of A along line L on 2: = 
and z = 1 are equal and opposite. 



A 2D mapping may be viewed as a discrete-time dy- 
namical system, and we will use mathematical methods 
developed for such systems. For more details see Guck- 
enheimer and Holmes— or Wiggins—. 

Analogous to an x-point in a 2D vector field, a fixed 
point xq of a 2D mapping is said to be hyperbolic when 
the eigenvalues As, A„ of the Jacobian matrix Jij = 
dFi^i/dFij at xq satisfy |As| < 1 < \Xu\- As with the x- 
point, the associated eigenvectors define linear subspaces, 
and the map Fi has global stable and unstable manifolds 
VF*(xo), IF"(xo) that are tangent to these linear sub- 
spaces at Xq. There are two branches of each manifold 
for each hyperbolic fixed point. 

By definition, M^*(xo) and M^"(xo) are invariant sub- 
spaces, meaning that if x S W^*(xo) then Fi(x) € 
^{xq), and similarly for W^{xo). Under the mapping 
Fi, points on W^^(xo) move closer to xo (along the curve), 
while those on W^"(xo) move further away. In the case 
of our magnetic field, W^*(xo) and 14^" (xq) correspond to 
curves on the boundary z = (or equivalently z = 1). 
Their invariance means that field lines starting on either 
manifold for z = must end on the same manifold for 
z = 1. The union of such field lines therefore defines 
a magnetic surface in the 3D domain generated by each 
manifold. 



A. Integrable fields 



IV. GENERAL FLUX PARTITION 

We propose a simple generalization of the 2D case: par- 
tition the flux in a 3D field by the hyperbolic fixed points 
of the field line mapping Fi{x,y) and their global man- 
ifolds. The physical nature of the partition is explored 
in this section, while the partition reconnection rate is 
defined in Section FVl 



The simplest type of 3D field to understand is an inte- 
grable field, where the field lines lie on a foliation of fiux 
surfaces. Figure[31[a) shows an example of such a field de- 
fined by adding a uniform z-component to a 2D magnetic 
field. Field lines of the 3D field lie on vertical surfaces 
that project on to field lines of the 2D field. The three 
null points of the 2D field now correspond to vertical field 
lines, thus to fixed points (ei, hi, h2). The scparatrices 
of the 2D field correspond to the global manifolds of the 



3D field line mapping, 
(a) integrable field 



(b) general field 




FIG. 3. Sketch of (a) an integrable field where field lines lie on 
flux surfaces, and (b) a more general 3D field created by per- 
turbing the integrable field. Thick black lines show magnetic 
field lines, including the three fixed points of the integrable 
field (hi, h2 hyperbolic and ei elliptic), which persist in the 
general field. Red and blue curves show where the global 
manifolds intersect the boundary. 



What is the flux of the "island" containing ei? There 
are two natural fluxes: (1) the vertical flux through the 
lower boundary within this region, and (2) the horizon- 
tal flux crossing the grey-shaded vertical surface between 
the ei and h2 field lines. The first flux does not exist in 
the original 2D field, but is measured straightforwardly 
from B • n on the lower boundary. The second flux is 
measured using the generalized flux function by the dif- 
ference A{ei) — A{h2). This is clearly analogous to the 
2D case (Section [ll]). Note that, in this integrable field, 
^(hi) = A{h.2), so the identical flux for this island re- 
gion would be measured by yl(ei)—^(hi). The "barrier" 
around the island comprises two magnetic surfaces: (1) 
a branch of Ty"(hi), which coincides with a branch of 
M^*(h2) (in red), and (2) a branch of VF*(hi), which co- 
incides with a branch of W^{h2) (in blue). 



ble manifolds. An intersection between a stable manifold 
and an unstable manifold of the same fixed point is called 
a homoclinic point, while an intersection between man- 
ifolds from different fixed points is a heteroclinic point 
(Figure |4]). In this paper, we shall not need to distin- 
guish between the two, and will refer to both as hetero- 
clinic points. 




FIG. 4. Example and notation for a heteroclinic tangle be- 
tween two hyperbolic fixed points hi and h2. The manifolds 
are curtailed to finite length for clarity. Pip q is a homo- 
clinic intersection. Pips pi, p2 are heteroclinic intersections, 
defining the lobe shaded in grey. 



The key result about heteroclinic intersections, first 
recognised by Poincare, is that a single intersection be- 
tween two manifolds 14^^ and W^ implies the existence 
of an infinite number of intersections between these same 
two curves. This simply follows from the fact that the 
intersection point lies on both manifolds. It cannot be a 
fixed point, and every iterate must also lie on both man- 
ifolds, by definition. The infinite number of intersections 
as one approaches either of the fixed points hi , h2 leads 
to a very convoluted path of the manifold curves. Called 
a homoclinic tangle, this is a major route to chaos in 2D 
mappings. This possibility of chaos in the field line map- 
ping is the major factor that complicates the partitioning 
of flux in a 3D flcld. 



B. Heteroclinic tangles 

Unfortunately, the simplicity of the integrable case be- 
lies the complexity typical of a general 3D magnetic field. 
If a small z-dependent perturbation is applied to Figure 
ini^a), the three fixed points will persist and maintain their 
elliptic/hyperbolic character, but the regular global man- 
ifolds will break down into heteroclinic tangles (Figure 
[3l^b)). In this generic situation, the stable and unsta- 
ble manifolds intersect trans versally at discrete points, 
rather than coinciding to form regular separatrices as in 
the 2D or integrable cases. It follows from the uniqueness 
of field lines that an intersection can take place only be- 
tween a stable manifold and an unstable manifold. Two 
stable manifolds can never intersect, nor can two unsta- 



C. Partial barriers 

Since the global manifolds for a 3D field can be in- 
finitely long (unlike in 2D), the regions of the fiux parti- 
tion must be defined by partial barriers: curves compris- 
ing segments of one or more global manifolds, ending at 
hyperbolic fixed points^^. 

To formally define a partial barrier, let W^[xi,X2] de- 
note the segment of W^ between two points xi , X2 . Con- 
sider an intersection point p G W^^ fl W^^ . This point 
p is a primary intersection point or pip if the segments 
W^ [hi,p] and W^^ [p, h2] intersect only at p (and pos- 
sibly at hi if h2 = hi; Rom-Kedar et al.—). A partial 
barrier starts and ends at hyperbolic fixed points (possi- 



bly the same) , and comprises one or more global manifold 
segments intersecting at pips. It includes no further fixed 
points. 

Figure [5ja) shows a partial barrier between hyperbolic 
fixed points hi and h2, with two segments W^ [hi, p] and 
W^ [p, h2] intersecting at pip p. The barrier separates 
the shaded region A from the unshaded region A'. There 
is nothing special about this choice of pip: choosing a 
different pip would redefine the barrier and also the shape 
of regions A and A'. But the partition fluxes are defined 
only by A at the fixed points, so are independent of the 
choice of partial barrier. 

(a) 




FIG. 5. A partial barrier and the turnstile lobes. 



The barrier in Figure [5] is called "partial" because cer- 
tain field lines cross it in the mapping Fi. While no 
magnetic field line may cross the magnetic surface gener- 
ated by each global manifold, this does not prevent field 
lines from crossing the partial barrier if it is made up 
of more than one global manifold. In the remainder of 
this section we show that the flux crossing the barrier in 
each direction under Fi is well-defined (independent of 
the choice of pip p) , and further that the net flux crossing 
the barrier is simply .4(h2) — .4(hi). 

The key to understanding which field lines cross a par- 
tial barrier is lobe dynamics^^'^—. A lobe is a closed re- 
gion bounded by the segments iy^Jpi,p2], W^hJPi:P2] 
between two adjacent pips pi, p2 (e.g., Figure]^. The 
important dynamical rules governing lobes are^^: 

I. Lobes map to lobes under Fi. This follows from 
continuity of the mapping and the fact that W^ and 
W are invariant manifolds that field lines cannot 
cross. 



2. Ordering of points on W" and W is maintained, 
so for a given pair of intersecting manifolds, there 
are a fixed number m of lobes lying between p and 
Fi(p), the same for any pip p. 

In our case, Fi is orientation-preserving (| J| > because 
-Bz > 0), so TO must be even. 

Consider again Figure [SJ where m = A. In the mapping 
Fi, the two lobes £'i, £'2 cross from Ato A' , while the two 
lobes Ci, C2 cross from A' to A. These four lobes, which 
are precisely those lying between p and Fi(p), are the 
turnstile lobes: they contain exactly those points which 
cross the partial barrier under Fi. The flux in a lobe 
L is measured by integrating $(£) = /, Bz{x, y, 0) dxdy 
on 2; = 0. What happens if we choose a different pip 
to define the partial barrier? The turnstile would then 
comprise different lobes. But the dynamical rules above 
guarantee that there would still be four turnstile lobes, 
and their fluxes would be the same as for the original 
choice of pip. 

Theorem 1 (Net flux). Let p be a pip of W^ and 

W^ defining a partial barrier between regions A and A' 
(oriented as in Figure\B^. For i ~ I,..,7Ti/2, let Ei be 
the lobes mapped from A to A' by Fi, and Ci the lobes 
mapped from A' to A. Then 

m/2 

^(h2) - ^(hi) = Y. (*(^^) - *(^'))' 



and this sum is independent of the choice of pip p. 




FIG. 6. Sketch of the magnetic surfaces in the 3D domain 
generated by field lines from W^i (in red) and Wy^^ (in blue). 
Some important magnetic field lines are shown in black. 



Proof. We already know from the rules of lobe dynamics 
that the sum is independent of the choice of pip. 

The sketch in Figure [6] illustrates the magnetic surfaces 
generated by W^ and W^^ , for the barrier in Figure [S] 
To derive our result, we consider two closed loops, one 



lying on each of these surfaces. Start with the following 
closed loop on the surface generated by W^ : 

Lu = F4p)UM/h"jFi(p),hi]UF7i(hi)UM/h"Jhi,p]. (7) 

Here the notation F^ (p) means the field line traced from 
p on the lower boundary to Fi (p) on the upper boundary, 
and F^^ means a field line traced downward from the 
upper boundary to the lower boundary. 

Now form a closed loop on the W^ surface: 

L, = F,(p)UTy4[Fi(p),h2]UF;i(h2)UM/^Jh2,p]. (8) 

The integral of A around each loop must vanish, so, using 
periodicity of A x e^. 



Aip)- [ A- dl-Aih 



Aip) 



'M'h2[P'Fi(P)l 

Subtracting (jlOp from ^ yields 



i)=0, (9) 

A- dl-A{h2)=Q. (10) 



^(h2)-^(hi) 



A-rfl- 



w. 



- [p,Fi(p)] 



A-dl. 



W^Jp,Fi(p)] 



(11) 

The right-hand side is the magnetic flux through a closed 
loop in the plane z = (or z = 1) encircling all of the 
turnstile lobes. Since the Fi{Ei) are encircled anticlock- 
wise and the Fi(Ci) are encircled clockwise, the result 
follows from Stokes' theorem. D 

Theorem [T] shows that, for two hyperbolic points con- 
nected by a partial barrier, the difference in their values 
of A is precisely the net flux crossing this partial barrier. 
If h2 = hi, i.e. the two manifolds belong to the same 
fixed point, then there can be no net flux across the bar- 
rier. In the limiting case that hi and h2 are connected 
by a regular separatrix (i.e., the two global manifolds 
coincide exactly, as in a 2D fleld), there are effectively 
infinitely many lobes of zero area. In this limiting case. 
Theorem [1] reduces to ^(h2) = ^(hi), as in 2D. 

For simplicity, our illustrations avoid secondary inter- 
sections between lobes (see Rom-Kedar et al.—). How- 
ever, even in the presence of secondary intersections, one 
may show that ^(h2) — -4(hi) gives the net flux across 
the partial barrier. 



D. A region of the general flux partition 

To get a feeling for the nature of our new flux partition, 
let us examine the region R in Figure[7lja), whose bound- 
ary is a chain of partial barriers comprising alternating 
segments of W^ and W from four hyperbolic fixed points 
hi, i = 1,...,4. The exact definition of each partial 
barrier is non-unique, owing to the freedom of choice of 
defining pip, but this does not affect any of the fluxes 
that we will describe here. 




FIG. 7. A region R of the general flux partition (a), and 
decomposition into subregions Ri (b). 



Since R is simply connected it must contain an elliptic 
fixed point e, because the topological degree of Fi on 
i? is 1, by definition of the boundary. In this respect, 
the situation is analogous to a similar region in a 2D 
magnetic field, which must contain an o-point (Figure[T|). 
However, there is a key difference when we try to define 
a "flux" of the region R. In the 2D field, the partial 
barriers would be replaced by regular separatrices, and 
the fiux function A would have the same value at each 
hi. The flux of the region would then be unambiguously 
defined as ^(e) — ^(hi). In the 3D case, such a unique 
flux cannot be defined. 

To see this, consider the differences in A between each 
pair of fixed points, which define 8 physical fluxes, shown 
in Figure [Tl^b) : 



^Pi = A{h2)-A{hi), 

^-2 = Aihs) - ^(h2), 

^3 = ^(h4) - ^(ha), 
i;4 = A{hi)-A{hi), 



h^A{e)-A{hi), 
h = A{e) - A{h2), 
h = A{e) - ^(ha), 
bi = ^(e) - ^(h4). 



The il^i are the net fluxes through each partial barrier, 
while the (j)i measure fluxes across surfaces in the domain. 
But these 8 fluxes are not all independent. Dividing R 
into 4 subregions Ri, as in Figure [TKb), the net flux into 
each Ri must vanish (for a periodic field), so we have the 
constraints 



V'i = 01 - 0! 



+ ll 



for t = 1, 



(12) 



Summing all of these equations leads to X)i^i ~ ^j 'Ex- 
pressing conservation of flux in the full region R. It is 
clear from (fT2)l that differences between the 4>i relate to 
net flux through the partial barriers. If all ipi = 0, as 
in a 2D field, all of the 0^ must be equal, giving us our 
uniquely defined flux. But if any of the ipi are non-zero, 
there is no meaningful single flux in R. 

Interestingly, we see that a change in A{e) adds the 
same amount to each (pi. So although the structure is 
not simple enough to define a unique flux in R, there 
is a unique reconnected flux. This emphasizes that it 



is the values of A at fixed points that define our flux 
partition, not the individual fluxes ipij 4'ii which are not 
independent. Note that this change in .4(e) docs not 
affect any of the V'i: so it represents a purely local non- 
ideal event within the region R. By contrast, changes in 
.4.(hi) affect more than one region of the partition, giving 
reconnection in the usual sense. 



V. MEASURING RECONNECTION 

Measuring the partition reconnection rate for our gen- 
eral flux partition is straightforward using the generalized 
flux function A. As in the 2D case, the partition fluxes 
are deflned entirely by the values of A at flxed points. 
In a general 3D field line mapping with chaotic regions, 
we cannot uniquely define the regions of the partition, 
because the definition of partial barriers is non-unique. 
But to measure reconnection we require not a partition 
of space into regions but a partition of flux. This we have 
shown to be well-defined. 

In the same way as the 2D case, we can define a global 
partition reconnection rate by summing over hyperbolic 
fixed points h^: 



A$T 



E 



dA{h.i 



dt 



(13) 



Faraday's law then implies that 

w X V X A = -V{(l) + V), 



dA 
'dt 



(17) 



where 0(x, t) is the electrostatic potential. Following a 
field line at velocity w, the rate of change of A is 



dA 
'dt 



dA 
'dT 



d 
di. 



A dl 



w X V X A-i- v(w • a; 



(w- A-^-t/.)! ' 



Fi{x,y) 

(x-.y) 



(18) 

dl (19) 

(20) 



using (|17p . At a fixed point, Fi(xo) = xq, and we now 
choose the gauge 4> so that (/> = w • A. This ensures that 
A becomes an ideal invariant whenever R = 0, and leads 
to 



9^ 

—-+W-VA 
dt 



Fi(xo) 



E- dl. 



(21) 



In other words, the rate of change of A following a fixed 
point corresponds to the integrated parallel electric field 
along the fixed point field line. 



How docs this partition reconnection rate relate to the 
definition of reconnection using integrated Eu7 In the 
latter definition, reconnection can occur anywhere with 
non-zero E^^: a global reconnection rate is defined^ by 
identifying distinct reconnection regions as local maxima 
of i?| I . These sites need not coincide with fixed point field 
lines, so the global reconnection rates from the two meth- 
ods may differ. The example in Section IVTl will demon- 
strate this. In fact, the change in y^ at a fixed point 
corresponds to the integral of i?|| along the fixed point 
field line itself, as we now show. 

Let Xq be a fixed point of the field line mapping Fi, 
either hyperbolic or elliptic. We consider a general non- 
ideal evolution with an Ohm's law of the form 



E + V X B = R, 



(14) 



where R is any non-ideal term (for example, in resistive- 
MHD, R = j/cr). In our magnetic field, where there are 
no closed field lines, we may write R = S/tp -|- u x B, 
where 



V(F,(x,2/))= / Rdl 

J (x.y) 



(15) 



is integrated along a field line, and u = B x (R— V'0)/i?^. 
Letting w = v — u, we may write Ohm's law (fT4)) as 



E + vi^ X B = V^}, 



(16) 



indicating that the field lines are frozen-in with a velocity 
w. which differs in gencralii from the plasma velocity v. 



VI. EXAMPLE 

To illustrate the ideas developed above, we present a 
particular example of a 3D magnetic field. The basic 
field (Section IVIA[) is given by an analytical expression, 
and is chosen to be generic in that the field line mapping 
contains both regular and chaotic regions, and the global 
manifolds form heteroclinic tangles. In Section fVIBI we 
perform a simple experiment where a growing toroidal 
fiux ring is added to the basic field. This demonstrates 
how our flux partition responds to localized 3D reconnec- 
tion. 



A. Structure of the basic field 

The basic fleld comprises six isolated magnetic fiux 
rings, superimposed on a uniform vertical fleld B^ ~ \. 
The ith flux ring is derived from a vector potential 

A. = a,fc.exp {. (^-^^f^iv-y^f _ (i_fl)!) e., 

(22) 
where the centre of the ring is at (x^, y,j, Zj), the parame- 
ter ki controls the flux, and the parameters a^ and U con- 
trol the radial and vertical extents respectively^^. Thus 



B = Vx L|e, + |e^ + ^A,j 



(23) 




>■ o.c 




FIG. 8. The 6-roll magnetic field. The left panel shows a 3D visualization (isosurfaces at B^ + By = 0.003 identify the toroidal 
flux rings), while the right panel shows a Poincare return map with the calculated global stable/unstable manifolds shown 
in blue/red respectively. The greyscale image on the lower boundary of the 3D visualization shows the "color map" used to 
identify fixed points. 



where we choose the parameter sets 

(x„i = l,...,6) = (l,0,-l, 1,0,-1), 

{y^,i^ 1,...,6) == (0.5,0.5,0.5,-0.5,-0.5,-0.5), 

(z„i = l,...,6)-(-20,-12,4,4,12,20), 

(fe„ i = 1, . . . , 6) = (1, -1, 1.1, -1, 1, -l)fco 

with fco = 0.08, all a, = 0.3^2 and k = 2. Notice that 
one of the rings has larger | fc^ | : this creates an asymmetry 
in the flux partition leading to a net flux across certain 
partial barriers. 

Figure [51[a) illustrates this magnetic field. Since both 
A and B are periodic, we can iterate the field line map- 
ping Fi{x,y) to produce a Poincare plot (Figure [51Jb)). 
This reveals that the field is structured into six regular 
elliptic regions — corresponding to the {x, y) locations of 
the six flux rings — separated by bands of chaotic field 
lines. 

The 2D greyscale image on the base of Figure ^a) 
shows a "color map"-S4 of the direction of Fi (z, y) — (x, y ) . 
We used this to identify fixed points as locations where 
all four colors meet, and to identify the Poincare index of 
each fixed point from the surrounding color sequence^^^. 
The precise locations (Table |I| were found by Newton- 
Raphson iteration using the color map as a first guess. 
There are six elliptic fixed points, at the centre of each 
regular region, and two hyperbolic fixed points between. 
Table |T] also shows the value of the generalized fiux func- 
tion A{x, y) at each fixed point, calculated by numerically 



integrating A along the appropriate magnetic field line 
from z = -24 to z = 24. 



TABLE I. Fixed points in the basic field ([23] 



Point Poincare index 



y A{x,y) 





1.0058 


0.5007 


0.1191 




-0.0021 


0.5083 


-0.1189 




-1.0003 


0.5003 


0.1314 




-1.0061 


-0.5010 


-0.1191 




0.0024 


-0.5081 


0.1190 




1.0007 


-0.5007 


-0.1194 




-0.4623 


-0.0457 


0.0007 




0.4663 


0.0410 






ei 

62 
63 

64 
65 
66 

hi 
h2 



The red and blue curves in Figure [SJb) show the global 
manifolds of the fixed points hi and h2. These have 
been "grown" numerically up to a finite length using the 
method of Krauskopf and Osinga— (see also England 
et al.—). As expected in a generic 3D mapping, the 
manifolds do not describe regular separatrices, but have 
degenerated into heteroclinic tangles. In addition to the 
manifolds of hi and h2 , we also show the corresponding 
manifolds emanating from the six hyperbolic points at 
infinity (effectively on the boundaries of this plot). 

An enlargement of the partial barrier between hi and 
\i2 is shown in Figure [5]J a). By identifying a pip p and 



computing Fi(p), we find that m = 2 for this example: 
i.e., one turnstile lobe crosses the partial barrier in each 
direction. It is apparent that the areas of these two lobes 
are unequal, i.e., there is a net flux across this partial 
barrier. Numerical integration shows that the two lobe 
fluxes are approximately 0.0008 and 0.0001, and indeed 
A{h2) - A{hi) = 0.0007, thus verifying Theorem [U for 
this example. Repeating the calculation for the other 
partial barriers, we find that those connecting hi with 
the boundary each have a net flux of 0.0007, while those 
connecting h2 with the boundary each have zero net flux. 
This is consistent with the values of A for the two hyper- 
bolic points (at infinity, ^ = 0). One can think of a net 
chaotic flux of 0.0007 encircling the left-hand hyperbolic 
point, crossing all four of its attached partial barriers. 
(These barriers must all have the same net flux since Fi 
is area preserving.) 



B. Effect of an isolated reconnection region 

To illustrate several key properties of our general flux 
partition, we consider the effect of adding a gradu- 
ally strengthening seventh flux ring to the basic field. 
This models the topological effect of a localised three- 
dimensional (non-null) diffusion region, as modelled by 
Hornig and Priest— and studied in the framework of 
general magnetic reconnection^S. 

Specifically, 



B = V X 



y 



X 



Ea. 



(24) 



where A^ are the same as the basic field for i = 1, . . . , 6, 
and the new ring has parameters zy = 28, 07 = 0.1, Z7 = 
1, and ky = O.Olt. The dependence of k-j on time t causes 
a gradual increase in the ring's azimuthal magnetic flux 
from c&y = at t = to $7 = y/TTa'j^kT[t) at time t. We 
shall illustrate how the reconnection associated with this 
new flux ring affects the partition fluxes for two differing 
locations (x7, 2/7). 

Case 1: At fixed point ha {xj = 0.4663, y^ == 0.0410). 
The resulting perturbation is shown in Figure [101 The 
fixed point remains at the same position, and for small 
t remains hyperbolic, though the structure of the global 
manifolds underlying the chaotic region is altered. At 
t ~ 2.5, there is a pitchfork bifurcation: the original fixed 
point becomes elliptic and a pair of new hyperbolic fixed 
points are formed. The fiux <1>7 of the new flux ring has 
become strong enough to perturb the field and create a 
new elliptic region. 

Note that, throughout this evolution, the change in 
^(h2) is exactly equal to the rate of increase of $7 (Fig- 
ure [TUJ^d)). AU of the new flux is counted by our parti- 
tion reconnection rate, because the fixed point field line 
passes through the centre of the reconnection region. In 
general, a reconnection site will not be aligned with the 
fixed point field lines in this way, so the full imposed fiux 



will not be measured by our partition reconnection rate. 
However, as suggested by the bifurcation in this exam- 
ple, if enough flux is reconnected then the structure of 
the underlying field will be modified, creating new fixed 
points that subsequently measure the new flux. 

Case 2: On a partial barrier [xj = j/7 = 0). Here 
the flux ring modifies field lines in the lobes of a partial 
barrier (Figure E]). Figures IHlJb) and (c) show that the 
lobes grow as the ring flux increases. This implies that 
more flux crosses the partial barrier in each direction. 
But the net flux across the barrier remains invariant be- 
cause neither fixed point field line passes through the re- 
connection region, so that the fixed point locations and 
A values cannot change. Notice that, even though the 
mapping near hi or h2 is unperturbed, the nearby lobes 
change significantly, because the global manifold passes 
through the reconnection region. Numerical computa- 
tion of the areas of the two turnstile lobes reveals that 
0.58(£'i -I- Ci) ss $7, roughly corresponding to a simple 
picture of the new flux being counted twice in the lobes: 
once in each direction across the partial barrier. But 
this reconnection has not changed the partition fluxes as 
defined by A at fixed points, hence the partition recon- 
nection rate is zero. 

To summarise, these examples illustrate three impor- 
tant properties of an isolated reconnection region: 

1. The position of a localised reconnection region 
within the background field determines its topo- 
logical effectiveness: i.e., the extent to which it 
changes the partition fiuxes, as measured by A at 
fixed points. Indeed, the fiuxes in the partition can 
change only if a fixed point field line passes through 
the reconnection region. 

2. If the reconnected flux becomes large enough, the 
reconnection region may perturb the original field 
sufficiently that new fixed points whose field lines 
pass through the reconnection region are created. 
Thus it will become visible to the fiux partition. 

3. If field lines from the lobes of a heteroclinic tan- 
gle pass through the reconnection region, then both 
the paths of the corresponding global manifolds and 
the lobe areas may change due to the reconnection. 
While this may alter the amount of flux in the turn- 
stile lobes, the net flux across the barrier cannot 
change if the associated fixed point field lines do 
not pass through the reconnection region. 



VII. CONCLUSION 

We have proposed a method to define and measure re- 
connection in a 3D magnetic field stretching between two 
boundaries. The flux is partitioned using the global man- 
ifolds of hyperbolic fixed points of the field line mapping 
between the boundaries. Individual fluxes in the parti- 
tion are defined as differences between the values of a 
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FIG. 9. Enlargement of the central partial barrier in Figure [Sj^b), showing the lobes in (a) the basic field, (b) Case 2 (Section 
IVTB| at i = 2, and (c) Case 3 at t = 6. 
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FIG. 10. Perturbation of the Poincare plot and global manifolds in Case 1, at (a) t — (basic field), (b) t = 2 and (c) t = 4. 
Panel (d) compares the reconnected flux measured using .4(h2) (thick solid line) to the imposed flux $7 (grey dashed line): the 
two curves coincide. The thin lines show the A values at the two new hyperbolic points following the bifurcation at f « 2.5. 



generalized flux function A{x,y) at fixed points. This 
is a natural generalization of the flux function in a 2D 
magnetic field, and maintains the key advantage that the 
reconnection rate (with respect to this partition) is mea- 
sured simply by the rate of change of A at fixed points. 
The associated partition reconnection rate is unique, and 
straightforward to compute: the main computational ef- 
fort required is identifying the fixed points in the field 
line mapping at successive times. 

Petrisor— has previously recognised that reconnection 
may be related to values of an action function (essen- 
tially A) on hyperbolic orbits, though that analysis was 
limited to non-twist area-preserving maps. Our results 
are more general, and have been derived in a more physi- 
cal way relating directly to the magnetic field itself. The 
interpretation of A as an action integral brings out a 
deep connection with the magnetic structure. Gary and 
Littlejohn— have shown that A{x, y) is the action in a 
variational formulation leading to the equations of the 
magnetic field lines. In other words, given A, and assum- 



ing displacements Sk with A • i5x = at the end-points, 
the Euler-Lagrange equations that extremise A are the 
equations of the field lines. 

Future work will consider in more detail how our par- 
tition reconnection rate compares with reconnection de- 
fined using the integrated parallel electric field, such as in 
numerical simulations (e.g. Pontin et al.—). There are 
certainly differences: we have seen in Section IVll how the 
visibility of a reconnection region (as defined with Eu) to 
our fiux partition depends on its location. We interpret 
this as a difference in the topological effectiveness of the 
reconnection, from the point of view of the flux partition. 
One way to refine the partition would be to integrate A 
over more than one iteration of Fi, measuring the values 
of A at the larger number of periodic points. 

The example in Section IVIBI also demonstrated that 
the flux in each direction across a partial barrier may 
be much larger than the net fiux, and furthermore may 
change under a non-ideal evolution even if the net fiux re- 
mains constant. Our existing flux partition is insensitive 
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to such changes. However, we note that it is possible to 
measure the flux of an individual lobe with the forinula 



^L)^ £ [A{Fl{p,))-A{Fl{p^))] 



(25) 



k— — c 



where P2 and pi are the defining pips of the lobe. 
This may be proved using a similar geometrical argu- 
ment to Theorem [H and is a straightforward generaliza- 
tion of a similar formula for lobe area in area-preserving 
maps2iiMi2^. Taking account of chaotic regions in this 
way is likely to be of particular importance in fields where 
heteroclinic tangles fill large areas of the plane, with mul- 
tiple intersections. In this case, the structure of the field 
is dominated by chaotic regions. While the flux parti- 
tion may still be defined in the same way using A values 
at the fixed points, the geometrical interpretation is less 
clear and requires further investigation. 

Finally, a limitation of this work at present is our re- 
striction to periodic fields where Bz{x, y,l) = Bz{x, y, 0). 
In future, we hope to relax this restriction thanks 
to recent theoretical progress in aperiodic dynamical 



system; 



.36 



The notion of hyperbolic trajectories, and 



their stable and unstable manifolds, may be extended to 
this more general setting. 
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